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Abstract. Despite recent major advances, the opacity re- 
mains a source of substantial uncertainty in the calcula- 
tion of solar models, and hence of solar oscillation fre- 
quencies. Hence it is of substantial interest to investigate 
the sensitivity of solar structure to changes in the opacity. 
Furthermore, we may hope from the precise helioseismic 
inferences of solar structure to obtain information about 
possible corrections to the opacities used in the model cal- 
culation. Here we carry out detailed calculations of the 
influence on solar models of changes in the opacity, in- 
cluding also evolutionary effects. We find that over the 
relevant range the response of the model is approximately 
linear in the opacity change, allowing the introduction of 
opacity kernels relating a general opacity change to the 
corresponding model changes. Changes in the convection 
zone can be characterized entirely by the change in the ini- 
tial composition and mixing length required to calibrate 
the model. 

Key words: Sun: evolution - Sun: interior - Sun: oscil- 
lations - Sun: convection zone - Sun: neutrinos - stars: 
interior opacity 



1. Introduction 

Accurate frequency measurements of thousands of modes 
of solar acoustic oscillation provide detailed information 
on the solar interior. In order to use these frequencies to 
derive internal structure and dynamics of the Sun, it is 
crucial to understand and limit the uncertainties in the 
computation of solar models and mode frequencies. Fur- 
thermore, it is of considerable interest to investigate the 
sensitivity of solar structure to changes in the input phys- 
ical parameters and properties of the solar interior. One 
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of the important physical properties in solar model cal- 
culations is the opacity, which is intimately linked to the 
properties of the solar oscillations through its effects on 
the mean structure of the Sun. Here we investigate the 
sensitivity of solar structure to local modifications to the 
opacity. 

Several authors have investigated the effect of opac- 
ity on the solar models and oscillation frequencies. In an 
early paper, Bahcall, Bahcall & Ulrich (1969) studied the 
sensitivity of the solar neutrino fluxes to localised changes 
in the opacity and equation of state and concluded that 
the neutrino capture rates are more sensitive to the equa- 
tion of state. The effects of artificial opacity modifica- 
tions on the structure of solar models and their frequen- 
cies were also examined by Christensen-Dalsgaard (1988). 
Constructing static models of the present Sun with an 
enhanced opacity near the base of the convection zone, 
he showed that the changes in structure and frequencies 
are approximately linear even to an opacity change of 
60%. The linearity of the response of the model to opac- 
ity changes was later confirmed by Christensen-Dalsgaard 
& Thompson (1991). In a detailed investigation, Korzen- 
nik & Ulrich (1989) attempted to improve the agreement 
between the theoretical and observed frequencies of oscil- 
lations by determining corrections to the opacity through 
inverse analysis. They found that the opacity inversion can 
only partially resolve the discrepancy. In a similar analy- 
sis, Saio (1992) obtained opacity corrections by fitting to 
low-degree frequency separations and helioseismically in- 
ferred sound speeds at a few points in the model; he found 
that much of the discrepancy between the Sun and the 
model could be removed by opacity changes of up to 50 %. 
More recently, Elliott (1995) investigated the helioseismic 
inverse problem as expressed in terms of corrections to the 
opacity. He derived kernels relating the opacity differences 
to the changes in frequencies, on the assumption that the 
change in luminosity could be ignored; he proceeded to 
carry out inversions for the opacity errors, based on the 
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observed solar oscillation frequencies, and neglecting pos- 
sible changes in composition associated with the change 
in opacity. He found that the differences between observed 
and computed oscillation frequencies could be accounted 
for by opacity changes of up to about 5 %. 

Here we carry out a detailed investigation of the sen- 
sitivity of solar structure to localised changes in the opac- 
ity, both for static and evolutionary solar models. The 
ultimate goal, taken up in a subsequent paper, is to de- 
termine the opacity corrections required to account for 
the helioseismically inferred properties of the solar inte- 
rior. In general the opacity k is a function of density p, 
temperature T and composition {X^}. However, it is evi- 
dent that information about the properties of the present 
Sun cannot constrain the opacity in such generality. Thus, 
for simplicity, we consider only opacity modifications that 
are functions 5\ogK{T) of temperature alone, log being 
the logarithm to base 10. If the opacity correction is suf- 
ficiently small that higher-order terms can be neglected, 
the response of any quantity F related to solar structure 
can be expressed in terms of a differential kernel Kp as 



^ = j Kf{T) 6\ogn{T) dlogT 



Here for simplicity we estimate the kernels from 
(<5F/F)(To) 



(1) 



(2) 



^^^^°^~/51og«(r)dlogT 
where 5F is the change corresponding to a suitably change 
5\ogK localised at T = Tq. Relations similar to Eq. (|l]) 
form the basis of inverse analysis (Gough 1985): if 5F is 
the difference between the observed and theoretical fre- 
quencies, Eq. (|l|) can in principle be inverted to deter- 
mine corrections to the opacity in the model (Korzennik 
& Ulrich 1989; Saio 1992; Elliott 1995). The kernels also 
provide a powerful visualization of the response in a given 
physical quantity of the solar model to a small pertur- 
bation in the input physics. As an example of this we 
evaluate kernels relating opacity changes to the structural 
differences in the solar model, the neutrino fluxes, and the 
small frequency separation between low-degree modes. 

2. Procedure 

To illustrate the sensitivity of the models to opacity 
changes we have computed extensive sets of comparatively 
simple models. These were based on the OPAL opacity 
tables (Iglesias, Rogers & Wilson 1992), the Eggleton, 
Faulkner & Flannery (1973) equation of state, and nu- 
clear reaction parameters from Parker (1986) and Bahcall 
& Ulrich (1988). Convection was treated with the mixing- 
length formulation of Bohm-Vitense (1958). The neutrino 
capture rates for the ^^Cl and ^^Ga experiments were ob- 
tained with the cross-sections given by Bahcall (1989). 
The heavy-element abundance Z was 0.01895. All models 
were calibrated to a fixed radius (6.9599 10^° cm) and lu- 
minosity (3.846 10^'^ ergs" '^) at the assumed solar age of 



4.6 10^ years by adjusting the mixing-length parameter ac 
and the composition. Diffusion and gravitational settling 
were ignored. Details of the computational procedure were 
described by Christensen-Dalsgaard (1982). 

It should be noted that our models are somewhat 
simplified, particularly in the choice of equation of state 
and the neglect of settling, compared with present state- 
of-the-art solar models (for a review see, for example, 
Christensen-Dalsgaard et al. 1996). Thus the values of 
the model parameters quoted, e.g., in Table 1 cannot be 
regarded as representative of the actual solar structure. 
However, we are here concerned with the differential ef- 
fects on the models of changes to the opacity; for this pur- 
pose, the present simplified physics is entirely adequate. 

We have considered two types of models: static models 
of the present Sun and proper evolution models, evolved 
from a chemically homogeneous initial Zero-Age Main- 
Sequence (ZAMS) model. In the static models, the hy- 
drogen profile X{q) where q = m/M (m being the mass 
interior to the given point and M is the total mass of the 
Sun), was obtained by scaling by a constant factor the 
profile Xr(g), obtained from a complete evolution model: 

X{q) = xXM , (3) 

where x ^^-s adjusted to fit the solar luminosity. In the 
evolution models the calibration to solar luminosity was 
achieved by adjusting the initial helium abundance Yq- 

The sensitivity of solar structure to opacity was inves- 
tigated by increasing k in a narrowly confined region near 
a specific temperature Tq according to 

log K = log + f{T^) , (4) 

where kq is the opacity as obtained from the opacity ta- 
bles; To was varied over the temperature range of the 
model. The function /(Tq) has the form 

2' 



/(To) = Aexp 



log T - log To 



(5) 



where the constants A and A determine the magnitude 
and width of the opacity modification. Henceforth the 
models computed with the same input physics but without 
opacity modifications will be referred to as reference mod- 
els whereas the perturbed ones will be referred as modified 
models. 

The detailed response of the oscillation frequencies to 
the opacity modifications will be considered separately. 
However, the importance of the small frequency separation 
as a diagnostics of the solar core motivates that we include 
it here, for comparison with the response of the neutrino 
flux. We characterize the separation, averaged over radial 
order n and degree /, by the parameter Dq defined by 

{Vnt - l^n-l,£+2)n,l ^ (4^ 4^ 6)Dq , (6) 

where asymptotically 



Do 



1 



dc dr 



47r2(no + 1/4 + /3) Jq dr r 



(7) 
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Table 1. Properties of static solar models. Here Xq is the initial hydrogen abundance, etc is the mixing-length 
parameter, dcz/R is the depth of the convection zone in units of photospheric radius R, Tc and Xc are central 
temperature and hydrogen abundance of the model of the present Sun, and the last three columns give the neutrino 
capture rates in the '^'^Cl and ""^Ga experiments, and the flux of neutrinos. 

Neutrino flux 



Model A Xo d^JR =*^C1 ^^Ga *B 
(lO^'K) (SNU) (SNU) lO^'cm'^s"^ 



SO 


0.0 


0.6978 


1.8098 


0.2739 


15.604 


0.3383 


8.092 


132.02 


5.747 


SA 


0.02 


0.6972 


1.8038 


0.2733 


15.601 


0.3380 


8.073 


131.97 


5.730 


SB 


0.1 


0.6948 


1.7813 


0.2709 


15.589 


0.3369 


8.006 


131.84 


5.670 


sc 


0.2 


0.6920 


1.7561 


0.2681 


15.576 


0.3355 


7.939 


131.72 


5.608 



Here uq is a suitable reference order, /3 is a constant pre- 
dominantly related to the structure of the outer parts 
of the model, and c is the sound speed. Details of the 
calculation of the oscillation frequencies and the evalua- 
tion of the average in Eq. (|^) were given by Christensen- 
Dalsgaard & Berthomieu (1991). The flt to the frequencies 
included modes with 1 = — 4 and radial order such that 
17 < n + 1/2 < 29; no was chosen to be 23. 

3. Results 

3.1. Linearity of the response 

The linearity of the response of model parameters to the 
opacity modiflcations was tested by constructing static 
models with A = 0.02, 0.1 and 0.2 (note that A = 0.1 cor- 
responds to a maximum change in opacity of 26%). The 
modiflcations were centred at logTp = 7, and the width 
was A = 0.02 (here and in the following T and To are 
measured in K). The properties of the reference and mod- 
ified models are presented in Table 1. It is evident that the 
changes in opacity required fairly substantial changes in 
the composition and mixing-length parameter to calibrate 
the model to obtain the correct luminosity and radius, and 
led to a significant change in the depth of the convection 
zone. 

Fig. 1 shows the change in sound speed, temperature 
and density resulting from the opacity change with A = 
0.1 and 0.2, the differences for A = 0.2 being scaled by 
1/2 to test linearity. Even though the curves do not co- 
incide precisely, the model changes are approximately lin- 
early related to even quite substantial changes in log k, 
as was also found in earlier investigations. We have in 
addition confirmed that the changes are relatively insensi- 
tive to the width A, provided the integrated modification 
/ <51ogKdlogT is kept fixed. 

A requirement for linearity is that the model does not 
become convectively unstable at the location where the 
opacity is modified. This constraint can approximately be 
formulated as 

{S In «;)„,ax = (In 10) A < ^ - 1 , (8) 
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Fig. 1. Differences at fixed fractional radius r/R between 
static modified models SB and SC and the reference model 
SO, in the sense (modified model) - (reference model). 
The heavy lines show results for A = 0.1 and the thin 
lines results for A = 0.2, multiplied by 1/2 to illustrate 
the extent to which the response is linear. The following 
differences are illustrated: 6r Inc^ (solid line); 6r hip (long- 
dashed line) ; and Sr In T (dash-dotted line) , In being the 
natural logarithm 

where V = dlnT/dlnp and Vad is its adiabatic value, 
both evaluated in the reference model. We note that 
this constraint is only barely satisfied for A = 0.2 and 
log To = 7: had we used A = 0.3 instead, the region around 
the modification would have become convectively unstable 
and the structure of the model would have been drastically 
altered. 

In the following we consider modifications computed 
using A = A = 0.02, corresponding to Model SA in Ta- 
ble 1; this is sufficiently localized to represent the linear 
response of opacity modifications confined to very nar- 
row local regions, similar to delta functions. Furthermore, 
this value of A ensures that the constraint (^ is satisfied, 
except in a few cases in the core during early stages of evo- 
lution, of little significance to the structure of the model 
of the present Sun. Specifically, a small short-lived ini- 
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tial convective core appeared in the evolution calculations 
for 7.11 < log To < 7.15; to avoid problems with non- 
linearities we simply suppressed the mixing of the compo- 
sition of the core in these cases. 



3.2. Model differences 

Changes in static solar models in response to opacity mod- 
ifications at various locations are illustrated in Fig. ^, 
in terms of differences (5^) between the modified and 
reference models at fixed fractional radius r/R. Corre- 
sponding results for evolution models are shown in Fig. ^. 
The changes in the initial hydrogen abundance and 
mixing-length parameter ac required to calibrate the mod- 
els are shown in Fig. || as functions of log Tq, together with 
other overall properties of the models. It is apparent that 
5r In T , Sr In c and 6r In p behave almost like step functions 
at the location T = Tq of the imposed opacity change; 
furthermore Sr In T and Sr In c are largely localized near 
this point. On the other hand, Srlnp has a gentler be- 
haviour. These properties are investigated in more detail 
in Sect. 



4.1 



below. 

Interestingly, the central temperature Tc depends in 
a non-monotonic fashion on logTp, SlnTc being even 
slightly negative for log To ~ 7.0 (see Fig. ^). The rea- 
son for this behaviour is not obvious, although it may 
be related to the steep decrease in Sr In T with increasing 
T, and hence decreasing r, below the region of immediate 
opacity modification (see also Sect. 4.1). The generation of 
the high-energy neutrinos is strongly coupled to Tc, and 
indeed the changes in the neutrino fluxes closely follow 
SlnTc (cf. Fig. I^b). In particular, there is also a slight 
decrease in the neutrino flux for an opacity increase near 
log To ~ 7. This is in contrast to the calculation of Bah- 
call et al. (1969) where the flux in models with increased 
opacity is higher than in the reference model at all values 
of To (cf. Fig. 1 of their paper). However, we note that 
Bahcall et al. used an opacity modification with a width 
in logT approximately six times as large as in our case; 
this is likely to have suppressed the finer structure in the 
dependence of the neutrino flux on Tq. 

Opacity increases generally lead to a decrease in the 
hydrogen abundance. Indeed, an increase in opacity would 
tend to decrease the luminosity; it follows from homology 
scaling (e.g. Kippenhahn & Weigert 1990) that this can 
be compensated for by an increase in the mean molecular 
weight, i.e., a decrease in Xq. When the opacity change 
is localized to the outer parts of the radiative region, the 
changes in the core of the model, and hence in the value 
of Xq required for calibration, become very small. 
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Fig. 2. Differences between modified static models and 
the reference model SO in the sense (modified model) - 
(reference model) . The modified models correspond to the 
following locations of the opacity increases: a log To = 7.0; 
b log To = 6.8; c log To = 6.7; d log To = 6.36. The vari- 
ables shown are: Srlnc^ (solid line); Sr\np (short-dashed 
line); Sr\np (long-dashed fine); (J^lnT (dash-dotted line); 
and SrX (dash-triple-dotted line) 
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Fig. 3. Differences between modified evolution models 
and the corresponding reference model, in the sense (mod- 
ified model) - (reference model). The modified models cor- 
respond to the following locations of the opacity increases: 
a log To = 7.0; b logTo = 6.7. The line styles have the 
same meaning as in Fig. ^ 
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Fig. 4. Responses of various global model quantities, as 
a function of the location Tq of the opacity modification, 
shown in the sense (modified model) - (reference model). 
Thin lines correspond to static models and heavy lines to 
evolution models. The thin vertical line corresponds to the 
temperature at the base of the convection zone in the ref- 
erence model, a Relative changes in central temperature 
Tc (dash-dotted line), depth of the convection zone (dot- 
ted line), and mixing-length parameter ac (dashed line), 
initial hydrogen abundance Xq (solid line) and central hy- 
drogen abundance Xc (dash-triple-dotted line) . (Note that 
for the static models, the relative changes in Xq and Xc 
are identical.) b Relative changes in the '^^Cl and ^^Ga 
capture rates (solid and dashed lines, respectively) and 
the flux of *B neutrinos (dot-dashed line) 
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Fig. and a comparison of Figs g and show sub- 
stantial differences between the static and the evolution 
models in the changes in X in the core of the model. 
In the static case, 5rX is determined by the scaling in 
Eq. (^), together with the change in mass fraction q at 
fixed r. The overall effect is that the change is somewhat 
smaller in the core than elsewhere. For the evolution mod- 
els, on the other hand, the change in the initial abundance 
is modified by the change in the rate of nuclear burning, 
which in turn depends predominantly on the change in 
temperature. As a result, the change in the central hy- 
drogen abundance Xc generally reflects the change in the 
central temperature Tc, as is evident from Fig. ^. This 
difference between static and evolution models in the be- 
haviour of the hydrogen abundance in the core leads to 
differences in the core response for other model variables, 
as is clear from a comparison of Figs || and |^; on the other 
hand, there is striking similarity between the results for 
static and evolution models elsewhere, indicating that the 
use of static models is indeed a useful and computation- 
ally efficient technique for investigating the response of 
the models to changes in the physics. 

The behaviour of Sr Inp and 6r Inp, particularly in the 
convection zone, shows a striking variation with Tq, the 
differences in the outer parts of the model changing sign 
for log To ~ 6.74. As discussed in Sect. 4.2 below, these 
changes are related to the complex behaviour of etc and 
result from the matching between the convection zone and 
the radiative interior, and the calibration of the model to 
fixed radius. 

We finally note from Fig. |4| that, as expected, the 
model differences are negligible for opacity modifications 
confined essentially to the convection zone, corresponding 
to log To ^ 6.25; this evidently follows from the fact that 
the structure of the nearly adiabatic part of the convection 
zone is independent of opacity. 



4. Analysis of the model response 

In general, the response of the model to changes in aspects 
of the physics is quite complex, as is also apparent from the 
results presented so far. However, it is possible to obtain 
some understanding of important aspects of the model 
changes from relatively simple analyses. Here we consider 
the response in the vicinity of the imposed opacity change, 
and the behaviour of changes in the convection zone. 



4.I. Local response of temperature and pressure 



The numerical results presented in Sect. 3^ indicate a 
great deal of regularity in the response near the point 
where the opacity modification is made. Some properties 
of this behaviour can be understood quite simply in terms 



of the perturbed equations of stellar structure. For sim- 
plicity, we assume that the change in local mass and lumi- 
nosity can be neglected; the latter assumption is certainly 
satisfied outside the nuclear-burning region, where the lu- 
minosity is fixed by the calibration, while we have found 
from the numerical results that the perturbation in mass 
is comparatively small. Also, we evidently consider only 
the radiative part of the star. Under these assumptions 
the model response is determined by just the equations 
of hydrostatic equilibrium and radiative energy transport 
which, upon linearization, yield 



dlnp 

dr 



d SrP 
dr p 
and 

d SrT _ dlnT 
dr 



T 



dr 



SrP 

P 
Sk 



6rP 
P 



(9) 



+ {Kp + 1) 



SrP 



int 



+ {kt - 4) 



SrT 

~T~ 



(10) 



In Eq. (|l0| ) Kp = {d\nK/d\n p)t and kt ~ (9 In k/9 In T)p, 
and {Sk/k)^^^ is the intrinsic opacity modification, corre- 
sponding to Eq. (|). In Eq. we neglected the effect on 
the opacity of the change in composition, required by the 
calibration of the model; this change is also found to be 
small, compared with the dominant effects. Still neglect- 
ing effects of composition and assuming the ideal gas law, 
the modifications in pressure, density and temperature are 
related by 

SrP ^Srp_^ SrT 

p p T ' 

It is evidently possible to carry out a complete descrip- 
tion of the modifications to the model through numeri- 
cal solution of these equations, with appropriate bound- 
ary conditions; the analysis should then include also the 
changes in mass and luminosity, as described by the lin- 
earized versions of the relevant equations, and the effect 
on the composition. Here, however, we approximate the 
equations further, in order to obtain a rough analytical 
solution which may be used to interpret the numerical 
results. From Eq. ( |lO| ) it is evident that an intrinsic mod- 
ification ((5K/K)int approximating a delta function induces 
a step function in SrT/T, while, from Eqs and (pi]), 
Srpjp results from an integration over this step function. 
Thus it is plausible, as confirmed by the numerical re- 
sults, that near the modification SrT jT is substantially 
larger than Srp/p- Consequently we neglect Srp/p in Eq. 
(|l|) to obtain Srp/p ^ -SrT/T and hence, from Eq. 



d SrT 

dlnT~ ~^ 
where C — 5 



.SrT 

T~ 



c- 



Sk 



(12) 

and kt are 

approximately constant, and take as boundary condition 



mt 

We assume that 
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that 5rT/T = {5rT/T)s at a suitable reference tempera- 
ture Ts <Tq. Then the solution to Eq. (|l|) is 



SrT 

where 



= < 



5rT 
SrT 



dlnT 



for T<To 



T 
for 



T > Tf 



(13) 







(14) 



int 



The exponent ( is large: 
lar radiative interior, kt 



at a typical point in the so- 
~ —2.5 and Kp ~ 0.5, so that 
C ~ 8. Consequently, the term in {Ts/T)i decreases very 
rapidly with increasing depth and, as a first approxima- 
tion, 6rT/T ~ KoiTo/T)< for T > Tq and zero otherwise. 
To estimate Srp/p we write Eq. (||) as, using Eq. (pT|), 

d SrP _ ^^l5rT 

T ' 



(15) 



dlnT p 

where V — dlnT/dlnp. Assuming V to be roughly con- 
stant over the relevant region and using the approximate 
solution for <5rT/T, we obtain 

' 5rP^ 



SrP 

P 



P 

SrP 

p 



CV 



T 



for T<Tn 



for T > Tn 



(16) 



where {Srp/p)o is the value of Srp/p aX T — Tq. 

To test these expressions, we consider in Fig. ^ the 
case, already presented in Fig. of a static model mod- 
ified with logTo = 6.8; since the analytical solutions de- 
pend predominantly on temperature, we show the results 
as a function of logT, rather than r/R. The analytical so- 
lution for SrT /T is clearly in good agreement with the nu- 
merical results; in particular, the rapid decrease, as T~'», 
for T > Tq accounts for the localized nature of the tem- 
perature modification. The fit is rather less satisfactory 
for Srp/p, although the analytical solution approximately 
recovers the magnitude of the variation in the vicinity of 
T = To; in fact, it is plausible that the neglect of the finite 
extent of (5Ac//«)int, the neglect of the change in composi- 
tion and the assumption that V is constant are of doubt- 
ful validity in this case. Even so, the analytical solutions 
do provide some insight into the nature of the model re- 
sponse. 

We finally note that the change in sound speed is given 

by 



Sr In c? — Sr In 7i 
~ Sr In 7i 



Sr In 



(5r In T — (5r In /i , (17) 

where in the last equality the ideal gas law was assumed; 
here 71 = {dhip/ d\a p)s, the derivative being taken at 
constant specific entropy s, and /.t is the mean molecular 
weight. Outside ionization zones of abundant elements. 
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Fig. 5. Comparisons of actual model differences (solid 
lines) with analytical approximations (dashed lines), for 
a static model with log To = 6.8, A = 0.02 and A = 0.02; 
to illustrate the behaviour of the solution, position in the 
model is indicated by logT. a Relative difference SrT jT 
in temperaturejthe analytical solution was obtained from 
Eqs (|l|) and (|l4|), neglecting the term in {SrT/T)^. b Rel- 
ative difference Srp/p in pressure; the analytical solution 
was obtained from Eqs ([l6|), fitting {Srp/p)o to the nu- 
merical solution at the location T = To of the opacity 
modification 

Sr In 71 can be approximately neglected; also, Sr In fi is al- 
most constant, except in the core where the composition 
has been altered by nuclear burning. As a result, as a first 
approximation Sr In is obtained from Sr In T by a con- 
stant shift, as is indeed observed in Figs || and ^. 



4-2. Model differences in the convection zone 

The bulk of the convection zone is very nearly adiabat- 
ically stratified; hence its structure depends only on the 
equation of state, the composition and the (constant) spe- 
cific entropy s. Indeed, it is evident from Fig. ^ that the 
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model is insensitive to the opacity in the convection zone. 
Assuming that the physics of the atmosphere is given, s is 
determined, within mixing-length theory, by the mixing- 
length parameter ac- If, as in the calculations presented 
in Sect. 3, we also assume that the equation of state and 
the heavy-element abundance Z are fixed, the structure 
of the adiabatic part of the convection zone is then es- 
sentially characterized by ac and the envelope hydrogen 
abundance Xg. Some simple properties of such convective 
envelopes are discussed in the Appendix. 

For model changes so small that the linear approxi- 
mation is valid, the difference in any quantity </> may be 
expressed in terms of the changes SXg and Sac m Xg and 



ac, as 



[dXg 



SXg 



dar 



Sac 



(18) 



The partial derivatives in Eq. (18) may be estimated from 
differences between envelope models differing only in Xg 
or ac as, for example {d(t)/dXs)a^-r — Sr4>/SXs. Here we 
consider three envelope models constructed using the same 
input physics as for the complete models described earlier, 
but with prescribed values of Xg and ac . Some properties 
of the models are listed in Table 2. In particular, K is the 
value of the co nsta nt in the adiabatic relation between p 
and p [cf. Eq. (A3) in the Appendix], averaged over the 
lower part of the convection zone, with r < 0.8R. The 
resulting derivatives, in the outer parts of the convection 
zone, are shown in Fig. ^ at greater depth the derivatives 
of In are very small, and the derivatives of In p are essen- 
tially constant, in accordance with the discussion in the 
Appendix [cf . in particular Eq. ( [A6| )] . 

For the models computed with changes in the internal 
opacity, the changes in ac and Xg arise from the calibra- 
tion to fixed radius and luminosity. More precisely, Xg 
is essentially determined by the condition that the lumi- 
nosity is fixed and, given Xg, ac is determined such as 
to ensure a continuous match at the base of the convec- 
tion zone, with a fixed surface radius. Given the partial 
derivatives determined from the envelope calculations, we 
may test Eq. (18) by comparing actual differences found in 
Sect. 3, resulting from opacity modifications, with those 
obtained from Eq. (|l^) with the changes in Xg and ac 
shown in Fig. Results for static models are shown in 
Figs I and |. It is evident from Fig. 0, for logTo = 7.0, 
that the simple relation (18) reproduces quite accurately 
the details of the sound-speed changes in the hydrogen and 
helium ionization zones. Also, the dominant contribution 
in this case clearly results from the change in ac- 

Fig. H shows the changes in squared sound speed and 
density at r = 0.99i?, as well as the contributions from the 
changes in ac and ATg, as a function of logTo. The depen- 
dence of the individual contributions clearly follows the 
variation of Sac and SXg, shown in Fig. In particular, 
the complex behaviour of Sac is refiected in a correspond- 
ing behaviour of Sr In and Sr In p, including the changes 




0.90 



0.92 



0.94 0.96 

r/R 



0.98 1.00 



s 



C 




0.90 



0.98 1.00 



Fig. 6. Estimated partial derivatives of model quantities 
in the outer parts of the convection zone, with respect to 
the envelope hydrogen abundance Xg (solid lines) and the 
mixing- length parameter ac (dashed lines), a Derivatives 
of squared sound speed, b Derivatives of density 

in sign for 6.7 < logTo < 6.8. For opacity changes closer to 
the base of the convection zone, the model differences are 
dominated by the change in ac, the contributions from the 
change in Xg becoming negligible. The main reasons for 
these variations are not obvious; however, it seems likely 
that the behaviour in the radiative interior (cf. Fig. |2|) 
causes the change of sign in of S^p/p and Srpj p at the 
base of the convection zone, as the location of the opacity 
modification moves closer to the surface. Since Srp/p is 
directly linked to the change in the polytropic constant K 
by Eq. (|A6|), and hence to the change in ac by Eq. ( |A10| ) 
this is the probable cause of the change in ac and hence, 



by Eq. (A12), in the depth of the convection zone (see also 
Fig. 1^). In turn, the depth of the convection zone largely 
controls the behaviour of Sr In T and Sr In in the outer 
parts of the radiative region: For logTo <i 6.74, the con- 
vection zone is shallower in the modified model than in 
the reference model. Consequently, the temperature and 
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Table 2. Properties of the envelope models. is the envelope hydrogen abundance, pcz is the pressure at the base 
of the convection zone and K is defined by the adiabatic relation (|A3| ) 



Model No. 








Pcz 


K 










( dyn cm^^) 


(cgs) 


Enl 


0.697867 


1.81586 


0.274105 


4.44746 10'^ 


9.52969 10"'' 


En2 


0.697867 


1.84273 


0.277000 


4.75692 10^3 


9.33427 10" 


En3 


0.702867 


1.81586 


0.273757 


4.34019 10" 


9.65974 10" 



0.0010 



0.0005 



S 0.0000 



-0.0005 




-0.0010 

0.95 0.96 0.97 0.98 0.99 

r/R 



1.00 



0.0002 




0.0000 



-0.0002 



6.2 



0.002 



0.000 



-0.002 



-0.004 




0.95 



0.96 



0.97 0.98 

r/R 



0.99 



1.00 



C 




6.2 



Fig. 7. Differences in squared sound speed (panel a) and 
density (panel b) in the upper part of the convection zone, 
for a static model modified at logTg = 7.0. The solid line 
shows the computed difference, from Fig. the dash- 
dot and dash-triple-dot lines show the contributions in 
Eq. (18) from Sa^ and SX^, respectively, and the dashed 
line shows the sum of these contributions 



Fig. 8. Differences in squared sound speed (panel a) and 
density (panel b) at r = 0.99i? in modified static models, 
as a function of the location log Tq of the opacity modifi- 
cation. The same line styles as in Fig. have been used 

interior which lies outside the opacity modification. This 
trend is reversed for logTo < 6.74. 



sound-speed gradients are less steep in the modified model 
just beneath the convection zone; since, as argued in the 
Appendix, the sound-speed difference is very small in the 
lower parts of the convection zone, this leads to a nega- 
tive sound-speed difference in that part of the radiative 



5. Tests of the opacity kernels 

As indicated by Eq. (1), the sensitivity of a physical 
quantity to the opacity perturbations can be measured 
by opacity kernels. In this section we present examples 
of such kernels. We furthermore test the accuracy with 
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logr 
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Table 3. Properties of solar models with reduced core 
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logr 



Fig. 9. Kernels Kc{r,T) relating the change in opacity 
to the sound-speed change (cf. Eq. |l^. Heavy lines show 
kernels based on evolution calculation, while the thin lines 
were obtained from static models, a Kernels as a func- 
tion of fractional radius r/R for logT = 6.4 (dotted line), 
logT = 6.6 (solid hue), logT = 6.8 (short-dashed hue), 
logT = 7.0 (dot-dashed line), and logT = 7.1 (long- 
dashed line). The upper abscissa shows logT, correspond- 
ing to the lower abscissa, b Kernels as a function of log T, 
for r/R = 0.70 (dotted line), r/R = 0.50 (solid line), 
r/R = 0.30 (short-dashed line), and r/R = 0.10 (long- 
dashed line). The upper abscissa shows r/R, correspond- 
ing to the lower abscissa 



which kernels can reproduce the results of a large opac- 
ity change, by applying them to solar models with ar- 
tificially reduced opacity values which simulate the ef- 
fect of Weakly Massive Interacting Particles (WlMPs) 
(Christensen-Dalsgaard 1992). 

Fig. 9 shows opacity kernels Kc{r,T) for sound speed, 
defined such that 



J Kcir, T)S log K{T)d\ogT 



(19) 



These kernels were determined from the model changes 
discussed in Sect. 3.2 and are illustrated as functions of 
r and of T. Results for both static and evolution models 
are shown. As discussed in Sect. 4, the kernel for a given 
temperature T is localized in r fairly close to the position 
in the model corresponding to that temperature. Also, it 
is interesting that the kernels can be determined quite 
accurately from static models, except in the core where 
the evolution of the hydrogen abundance is substantially 
affected by the opacity modification. This provides some 
justification for the procedure adopted by Elliott (1995), 
which neglected possible evolutionary effects. 

To illustrate the effects on important global properties 
of the model Fig. 10 shows kernels for the neutrino fluxes 
^i, as well as the small frequency separation GDq (cf. Eq. 
6). Due to the high sensitivity of the neutrino fluxes to 
the opacity modifications, the neutrino kernels are defined 
specifically by 

Sln^^ = J K^{T)SlogK{T)d\ogT . (20) 

It is evident that these kernels correspond precisely, 
apart from a scaling, to the neutrino-flux differences shown 
in Fig. 4b. Also, it is striking that even for the neutrino 
flux there is little difference between the results for static 
and evolution models. In contrast, kernels for the small fre- 
quency separation, which depends strongly on the compo- 
sition gradient in the core, cannot be estimated accurately 
from the static models. 

To test these kernels, we use the models of Christensen- 
Dalsgaard (1992), with reduced core opacity. Specifically, 
the opacity modification was determined from 

logK = log«o - A/ (logT), (21) 
where 

2' 



/(logT) = 



, lOgT-lOgTl , .rrr^rr 

[ AlogT jj'^f^<^^ 
otherwise; 

(22) 

as did Christensen-Dalsgaard (1992), we used AlogT — 
0.04 and logTi = 7.1. This form of / provides an opacity 
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Fig. 10. Kernels for global model quantities. Heavy lines 
show kernels based on evolution calculation, while the thin 
lines were obtained from static models, a Kernels for loga- 
rithmic changes in the ^^Cl and ^^Ga capture rates (solid 
and dashed lines, respectively) and in the flux of '^B neu- 
trinos (dot-dashed line) (cf. Eq. 20). b Kernels for the 
averaged small separation 6-Do (cf. Eqs |^ and |^) 

decrease over a well defined region in log T with a contin- 
uous transition to zero at lower temperature. (Note that 
in these models the unmodified opacity kq was obtained 
using the Los Alamos Opacity Library of Huebner et al. 
1977.) It was found that the neutrino flux was reduced by 
a factor of more than three for A = 0.4 (see also Table 3). 

Table 3 compares the values of the ^^Cl neutrino flux 
and the small frequency separation computed for the mod- 
els with reduced opacity with those estimated by means 
of the kernels shown in Fig. 10. It is evident that there is 
good agreement between the values of Dg for both A = 0.2 
and 0.4. The agreement in neutrino flux is somewhat bet- 
ter for A = 0.2, the error being 10 - 15 % for A = 0.4. 
However, it is remarkable that linearized expressions of 



Fig. 11. Differences at fixed fractional radius r/i? between 
models with opacity modification to simulate the effect 
of WIMPs (cf. Christensen-Dalsgaard 1992) and the cor- 
responding reference models, in the sense (modified) - 
(reference). The thick lines represent the values obtained 
from the kernels while the thin lines show the original 
differences. Variables shown are Srhic (solid line); Jrlnp 
(dashed line); 6r\nT (dash-dot line). Some properties of 
the models are provided in Table 3. a Differences between 
the modified static model SW2 and the reference model 
EWO. b Differences between the modified evolution model 
EW2 and the reference model EWO 

the form given in Eq. (^ remain reasonably precise for a 
reduction in opacity by more than a factor of two. 

To test the precision with which the kernels repro- 
duce the structure of the model. Fig. 11 compares differ- 
ences in and p obtained from the kernels in Fig. 9 and 
the corresponding density kernels, with the actual model 
differences. Results are shown both for the static model 
SW2 and the evolution model EW2, both corresponding 
to yl = 0.4. As before, the agreement is quite close both 
for evolution and for static models. 
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6. Conclusion 

Our understanding of the solar internal structure has im- 
proved significantly over the last few years due to the rapid 
progress in helioseismology as well as in the construction of 
solar models with improved input physics. Although much 
of the knowledge has been obtained through inverse anal- 
ysis, the starting point is always solar models; thus it is 
important to investigate the sensitivity of the solar struc- 
ture to the input physics and to determine the presence of 
any region which is particularly sensitive to a specific pa- 
rameter. Since the major uncertainties in the microphysics 
come from the opacity, we have focussed on the examina- 
tion of the solar structure by means of a localised opacity 
change as a function of temperature. 

The sensitivity of the solar structure was represented 
by kernels relating the opacity changes to neutrino flux, 
frequency separation at low degree as parametrised by Dq 
and the structure difference between modified and refer- 
ence models. These kernels were subsequently used to de- 
rive the same parameters corresponding to a reduction 
in opacity by a factor of more than two in the models 
of Christensen-Dalsgaard (1992), simulating the effects 
of Weakly Interacting Massive Particles. We found that 
the kernels were remarkably successful in estimating the 
changes in the solar structure caused by even such a large 
change in the input physics. 

A natural next step in this investigation is to study sys- 
tematically how the oscillation frequencies are affected by 
opacity changes and also how well these frequency changes 
can be reproduced by kernels. Furthermore, in a prelim- 
inary analysis Tripathy, Basu & Christensen-Dalsgaard 
(1997) found that much of the current discrepancy be- 
tween the helioseismically inferred solar sound speed and 
the sound speed of a standard solar model can be un- 
derstood in terms of modest modifications to the opacity. 
More detailed analyses of this nature are now under way. 
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Appendix A: Some properties of convective 
envelopes 

For convenience, we present some simple approximate 
properties of convective envelopes, which are useful for 
the interpretation of the numerical results in Sect. 4.2. 
Further details were provided by Christensen-Dalsgaard 
(1997). Note also that Baturin & Ayukov (1995, 1996) 



carried out a careful analysis of the properties of convec- 
tive envelopes and their match to the radiative interior. 

In the convection zone the stratification departs sub- 
stantially from being adiabatic only in a very thin region 
just below the photosphere. Thus the structure of the bulk 
of the convection zone is determined by the equation of 
hydrostatic equilibrium 
dp Gmp 
dr 



(Al) 



(A2) 



(G being the gravitational constant and m the mass inte- 
rior to r), and the relation 
1 _ dlnp _ 1 
7 d In p 7i 
To this approximation, the convection-zone structure 
therefore depends only on the equation of state, the com- 
position and the constant specific entropy s, the latter be- 
ing fixed by the mixing-length parameter Uc- In practice, a 
more useful characterization of the bulk of the convection 
zone follows from noting that 71 is approximately con- 
stant outside the dominant ionization zones of hydrogen 
and helium; then Eq. (A2) shows that p and p are related 

by 

p ^ Kp-^' , (A3) 

where the constant K is closely related to s. Assuming 
71 again to be constant and neglecting the ma ss i n the 
conv ection zone, it is readily shown from Eqs (Al) and 
( [A^ ) that u = pI p — c? is given by 

(A4) 

approximately valid in the lower parts of the convec- 
tion zone (e.g. Christensen-Dalsgaard 1986; Dziembowski, 
Pamyatnykh & Sienkiewicz 1992; Christensen-Dalsgaard 
& Diippen 1992; Baturin & Mironova 1995); here M is 
the mass of the model and R* ~ R, its surface radius. 



P 



1-1/71 ^J^-l/7l2:i \qm 



(A5) 



From Eqs (A3) and (A4) wc also obtain 

1 1 

71 \r R^ 

Eq. ( |A4| ) indicates that u, and therefore c, depend little 
on the details of the physics of the model; in particular, 
if M and R are assumed to be fixed, as in the case of 
calibrated solar models, 5rU ~ 0. This is confirmed by 
Figs ^ and H, which show that 5r In c is very small in the 
bulk of the convection zone. From Eq. (A5) it furthermore 
follows that 

5r Inp ~ (5r In p ~ 5 In K (A6) 

7i - 1 

are approximately constant, again in accordance with 
Figs 1^ and ^. Finally, assuming the ideal gas law, we have 
that 

5r hi r ~ In . ( A7) 

These relations may be used to investigate the changes 
in the model resulting from changes in the convection- 
zone parameters. Of particular interest are conditions at 
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the base of the convective envelope, defining the match to 
the radiative interior. Neglecting possible convective over- 
shoot, Vad = Vrad, the radiative temperature gradient, 
at this point. For calibrated models the luminosity L is 
unchanged; then Vrad Kp/T'^. Assuming also that Vad 
is constant and that the heavy-element abundance Z is 
fixed, we find that the change in the pressure Pcz at the 
base of the convection zone is related to the changes in K 
and the envelope hydrogen abundance Xg by 

6 In Pcz — — — z — — TT— 7z 7T SlnK 



(4- kt)(7i - 1) -7i(^p + 1) 
71 [(4 - fiT)px - iix] 



SlnXs (A8) 



(4- kt)(7i - 1) -7i(^p + 1) 
(see also Christensen-Dalsgaard 1997); here 

kp — {dhiK/dhip)T,x 7 kt = {dlnK/dhiT)px , 

Kx = {d\nK/d\nX)p^T , f^x = {d\nfi/d\nX)z ■ 

Having obtained Slnpcz, the change in the depth of the 
convection zone can be determined from Eq. (A5) as 



(Slndcz =i — ^[^Inif + (71 - l),51npez] 
7i 



(A9) 



where is the radius at the base of the convection 
zone. [Note that this relation differs from Eq. (11) of 
Christensen-Dalsgaard (1997). There it was assumed that 
the interior properties of the model were largely un- 
changed while the surface radius was allowed to change; 
this is the case relevant, for example, to the calibration of 
solar models to have a specific radius. Here we have kept 
R fixed and the change in dcz corresponds to changes in 
the properties of the radiative interior of the model.] 

It is of some interest to compare these simple relations 
with the numerical results obtained for the envelope mod- 
els listed in Table 2. To do so, we first need to relate the 
change in K to the changes in ac and Xg. From the results 
in Table 2 

As discussed by Christensen-Dalsgaard (1997), the rela- 
tion between ac and K, at fixed composition, follows sim- 
ply from the properties of the mixing-length theory; the 
result of such an analysis, using the properties of the ref- 
erence envelope model Enl, agrees quite closely with the 
value given above. It is less straightforward to derive a sim- 
ple expression for the relation between 6\nK and SlnXs- 
To evaluate the changes in pcz and dcz, we need the 
derivatives of k and fi. At the base of the convection zone 
in the reference envelope Enl we have the following values: 
kp = 0.63, kr = —3.70, kx = —0.16, and fix — —0.54. 
Thus, from Eq. (^) we obtain Slnpcz -3.20(5 In + 



2.996lnXs, and hence, from Eq. (|/ 
Slndcz ~ -0A9S\nK + 0.87S\nX, 



Using also Eqs (AlO) we find 
Slndcz ^ 0.69(5 In - 0.07(5 In AT, 



(All) 



(A12) 



For comparison, the results in Table 2 give 
'dlndcz\ / dlndcz 



din Uc 



0.72 



dlnXs 



0.18. (A13) 



Evidently, the Q!c-derivative is in reasonable agreement 
with Eq. (A12); although the agreement appears less sat- 
isfactory for the derivative with respect to X g, it should 
be noticed that small coefficient in Eq. (A12) arises from 
near cancellation between the contributions from the two 



terms in Eq. ( All ) 
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